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Abstract 

We perform an asymptotic analysis of general particle systems arising in collective 
behavior in the limit of large self-propulsion and friction forces. These asymptotics impose 
a fixed speed in the limit, and thus a reduction of the dynamics to a sphere in the velocity 
variables. The limit models are obtained by averaging with respect to the fast dynamics. 
We can include all typical effects in the applications: short-range repulsion, long-range 
attraction, and alignment. For instance, we can rigorously show that the Cucker-Smale 
model is reduced to the Vicsek model without noise in this asymptotic limit. Finally, a 
formal expansion based on the reduced dynamics allows us to treat the case of diffusion. 
This technique follows closely the gyroaverage method used when studying the magnetic 
confinement of charged particles. The main new mathematical difficulty is to deal with 
measure solutions in this expansion procedure. 
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1 Introduction 



This paper is devoted to continuum models for the dynamics of systems involving living 
organisms such as flocks of birds, school of fish, swarms of insects, myxobacteria... The 
individuals of these groups are able to organize in the absence of a leader, even when starting 
from disordered configurations ^37j. Several minimal models describing such self-organizing 
phenomenon have been derived |38l 1281 [TO] . Most of these models include three basic effects: 
short-range repulsion, long-range attraction, and reorientation or alignment, in various ways, 
see |33) and particular applications to birds [32] and fish [U [2] . 

We first focus on populations of individuals driven by self-propelling forces and pairwise 
attractive and repulsive interaction [34:\ [25] . We consider self-propelled particles with Rayleigh 
friction [17[ [T6] [TT| 114] . leading to the Vlasov equation in d = 2, 3 dimensions: 

dtf + V ■ V^f + a=(t, x) ■ V^,f + ^dw^{f{a - I3\vf)v} = 0, {t, x, v) €R+xR'^x (1) 

where = f^{t, x,v) > represents the particle density in the phase space (x, t;) G M'^ x 
at any time t G M+, a"^ stands for the acceleration 

a'{t,-) = -VccU*p'it,-), p'{t,-)= [ r{t,;v)dv, 

and U is the pairwise interaction potential modelling the repelling and attractive effects. 
Here, the propulsion and friction forces coefficients = j > 0, f3'^ = j > are scaled in 
such a way that for e — )• particles will tend to move with asymptotic speed -^J^- These 
models have been shown to produce complicated dynamics and patterns such as mills, double 
mills, flocks and clumps, see [25]. Assuming that all individuals move with constant speed 
also leads to spatial aggregation, patterns, and collective motion [21} [26], 

Another source of models arises from introducing alignment at the modelling stage. A 
popular choice in the last years to include this effect is the Cucker-Smale reorientation pro- 
cedure |2U] . Each individual in the group adjust their relative velocity by averaging with 
all the others. This velocity averaging is weighted in such a way that closer individuals in 
space have more influence than further ones. The continuum kinetic version of them leads to 
Vlasov-like models of the form ([T]) in which the acceleration is of the form 

a%t,-) = -H*rit,.), 

where * stands for the (x, t;)-convolution, abusing a bit on the notation, with the nonneg- 
ative interaction kernel H : M?'^ — > M^. In the original Cucker-Smale work, the interac- 
tion is modelled by H{x, v) = h{x)v, with the weight function h being a decreasing radial 
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nonnegative function. We refer to the extensive literature in this model for further details 

[311 [291 nana Eg. 

In this work, we will consider the Vlasov equation ([T]) where the acceleration includes the 
three basic effects discussed above, and then takes the form: 

a'it, •) = -V,U * p%t, ■)-H* fit, •) . (2) 

We will assume that the interaction potential U € C^(M'^), U bounded continuous with 
bounded continuous derivatives up to second order, and H{x,v) = h{x)v with h G C^iW^) 
and nonnegative. Under these assumptions the model (HI)-© can be rigorously derived as a 
mean- field limit |36l [HI [Ml [IDl [3] from the particle systems introduced in [251 [20] . 

We will first study in detail the linear problem, assuming that the acceleration a = a{t, x) 
is a given global-in-time bounded smooth field. We investigate the regime e \ 0, that is 
the case when the propulsion and friction forces dominate the potential interaction between 
particles. At least formally we have 

r = / + e/«+e¥') + ... (3) 

where 

div,{/(a-/3|7;p)t;} = (4) 

dtf + diY^v) + div,(/a(t, x)) + div,{/(^)(a - P\vf)v} = , (5) 

up to first order. Therefore, to characterize the zeroth order term in the expansion we need 
naturally to work with solutions whose support lies on the sphere of radius r := \J aj ji 
denoted by rS with S = {f G M'^ : 1^1 = 1}. In turn, we need to work with measure solutions 
to Q which makes natural to set as functional space the set of nonnegative bounded Radon 
measures on M*^ x M*^ denoted by x W^). We wiU be looking at solutions to ([1]) 

which are typically continuous curves in the space x M*^) with a suitable notion 

of continuity to be discussed later on. We will denote by /^(t, x, d(x, w) the integration 
against the measure solution f^{t,x,v) of ([T|) at time t. For the sake of clarity, this is done 
independently of being the measure f^(t) absolutely continuous with respect to Lebesgue or 
not, i.e., having a L^{W^ x M.'^) density or not. 

Proposition 1.1 Assume that (1 -1- |f G A^^(IR'^). Then F is a solution to 1^ if and 
only if supp F C {0} U r§. 

The condition ^ appears as a constraint, satisfied at any time t G M+. The time evolution 
of the dominant term / in the Ansatz ([3|) will come by eliminating the multiplier f^^^ in ([5|), 
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provided that / verifies the constraint In other words we are allowed to use those test 
functions 'il>{x,v) which remove the contribution of the term div„{/(^^(a — /3|t;p)t;} i.e., 



{a - P\v\^)v ■ V^ip P'>{t, X, v) d{x, v) = 0. 

<R<^ 

Therefore we need to investigate the invariants of the field (a — f3\v\'^)v ■ V„. The admissible 
test functions are mainly those depending on x and v/\v\,v 7^ 0. The characteristic flow 
(s, v) — )• V(s; v) associated to ^{a — (3\v\'^)v ■ V„ 

^ = \a-(3\Vis;vt)Vis;v), ViO;v)=v 
as e 

will play a crucial role in our study. It will be analyzed in detail in Section [3j Notice that the 
elements of {0} U rS are the equilibria of (a — /3|f p)u • V^,. It is easily seen that the jacobian 
of this field 

^^,{{a - P\v\'^)v} = (a - /3|t;|^)/ -2^v0v 

is negative on rS, saying that rS are stable equilibria. The point is unstable, dy{{a — 
P)f}|^=o = al. When e \ the solutions {f^)e concentrate on M.'^ x ({0}UrS), leading to 
a limit curve of measures even if {f'^)s were smooth solutions. We can characterize the limit 
curve as solution of certain PDE whenever our initial measure does not charge the unstable 
point 0. 

Theorem 1.1 Assume that a £ W'^'°°(R'^)), ° G X+(M'^xM'='), supp f C 

{{x,v) : \v\ > vq > 0}. Then {f^)e converges weakly -k in L°°(M+; 7Wf,(M^ x M'^)) towards the 
solution of the problem 

5J + div,(/t;)+div,|/(^/-^^a| =0 (6) 

div,{/(a-/3|^;|2)t;} = (7) 
with initial data /(O) = (/™) defined by 

i^{x,v) (f'') {x,v)d{x,v) = / il^ ( x,r-^] f''{x,v)d{x,v) , 



\V\ 



for all i; G C^{M.'^ x 

In the rest, we will refer to (/'") as the projected measure on the sphere of radius r 
corresponding to Let us point out that the previous result can be equivalently written 
in spherical coordinates by saying that f{t,x,uj) is the measure solution to the evolution 
equation on {x,uj) G M'^ x rS given by 

dtf + div,(/a;) + div^ // ( / _ ^ ) a I = . 



These results for the hnear problem, when a{t, x, v) is given, can be generalized to the non- 
linear counterparts where a{t,x) is given by ([2]). The main result of this work is (see Section 
[5] for the definition of Vi): 

Theorem 1.2 Assume that U G Cl{W^), H{ x,v) = h{x)v with h G C^(M'^) nonnegative, 
/in g ^^(j^d X ^d^^ g^pp Jin ^ {(x,v) : \x\ < Lo,ro < |t;| < Ro} with < tq < r < Rq < oo. 
Then for all 6 > 0, the sequence {f^)e converges in C{[6,oo);Vi(M.'^x'K'^)) towards the measure 
solution f{t,x,uj) on {x,u}) G M'^ x rS of the problem 

dtf + div,{f^) - div^ !^f (^I - ^{cj oj)^ {V^U*p + H*f)^ = (8) 

with initial data /(O) = (/™)- Moreover, if the initial data /™ is already compactly supported 
on Bl^ X rS, then the convergence holds in C(M-|_ ; Pi (M'^ x M^)). 

Let us mention that the evolution problem (El) on X rS was also proposed in the 
literature as the continuum version [23j of the Vicsek model [6S\ [TE\ without diffusion for the 
particular choice U = and H{x,v) = h{x)v with h{x) some local averaging kernel. The 
original model in [38\ [T8] also includes noise at the particle level and was derived as the mean 
filed limit of some stochastic particle systems in [3]. In fact, previous particle systems have 
also been studied with noise in [3j for the mean-field limit, in [30] for studying some properties 
of the Cucker-Smale model with noise, and in [22^ [27] for analyzing the phase transition in 
the Vicsek model. 

In the case of noise, getting accurate control on the particle paths of the solutions is a 
complicated issue and thus, we are not able to show the corresponding rigorous results to 
Theorems 11.11 and 11.21 Nevertheless, we will present a simplified formalism, which allows 
us to handle more complicated problems to formally get the expected limit equations. This 
approach was borrowed from the framework of the magnetic confinement, where leading order 
charged particle densities have to be computed after smoothing out the fluctuations which 
correspond to the fast motion of particles around the magnetic lines O [6l [71 [8] . We apply 
this method to the following (linear or nonlinear) problem 

dtf + dW.if'v} + divAfa} + ^div,{f{a - (3\v\^)v} = A,f (9) 

with initial data /^(O) = /'"^ where the acceleration a G L°°{R+;W'^'°°{R'^)) and /"^ G 
A4'^{M.'^ X M*^). By applying the projection operator (•) to ([H), we will show that the limiting 
equation for the evolution of f{t,x,oo) on {x,uj) G M'^ x rS is given by 

dtf + div,(/a;) + div^ I f f I - ^ oj)) a] = A^f (10) 



where IS.^ is the Laplace-Beltrami operator on rS. 

Our paper is organized as follows. In Section [2] we investigate the stability of the char- 
acteristic flows associated to the perturbed fields v ■ V^. + a ■ + ^{a — P\v\'^)v ■ V^- The 
first limit result for the linear problem (cf. Theorem II. ip is derived rigorously in Section [3l 
Section 2] is devoted to the proof of the main Theorem 11.21 The new formalism to deal with 
the treatment of diffusion models is presented in Section O The computations to show that 
these models correspond to the Vicsek models, written in spherical coordinates, are presented 
in the Appendix Rl 



2 Measure solutions 

2.1 Preliminaries on mass transportation metrics and notations 

We recall some notations and result about mass transportation distances that we will use in 
the sequel. For more details the reader can refer to [39| [T5]. 

We denote by Vi{R'^) the space of probability measures on M with finite first moment. 
We introduce the so-called Monge-Kantorovich- Rubinstein distance in ViiW^) defined by 



Wi{f,g) = sup 



^{u){f{u) - g{u))du 



.ip£ Lip (M'^), Lip ((^) < 1 



where Lip(M^) denotes the set of Lipschitz functions on and hip{ip) the Lipschitz constant 
of a function ip. Denoting by A the set of transference plans between the measures / and 
g, i.e., probability measures in the product space x M*^ with first and second marginals / 
and g respectively 

f{y)= / 7r(y,z)dz, g{z) = / TT{y,z)dy 

jRd JM.d. 

then we have 

Wi{f,g)=mi \ \y - z\7r{y,z)d{y,z)\ 

by Kantorovich duality. Vi{M.'^) endowed with this distance is a complete metric space. Its 
properties are summarized below, see|39j. 



Proposition 2.1 The following properties of the distance Wi hold: 

1) Optimal transference plan: The infimum in the definition of the distance Wi is 
achieved. Any joint probability measure tTq satisfying: 

Wi{f,g) = / \y - z\d7ro{y,z) 
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is called an optimal transference plan and it is generically non unique for the Wi- 
distance. 

2) Convergence of measures: Given {fk}k>i o.nd f in ViiW^), the following two as- 
sertions are equivalent: 

a) Wi{fk,f) tends to as k goes to infinity. 

h) fk tends to f weakly -k as measures as k goes to infinity and 



sup / \v\fk{v)dv^O as i? — )• +cxd. 

fe>l J\v\>R 



Let us point out that if the sequence of measures is supported on a common compact set, 
then the convergence in T^i-sense is equivalent to standard weak-* convergence for bounded 
Radon measures. 

Finally, let us remark that all the models considered in this paper preserve the total mass. 
After normalization we can consider only solutions with total mass 1 and therefore use the 
Monge-Kantorovich- Rubinstein distance in Vi{M.'^ x W^). From now on we assume that the 
initial conditions has total mass 1. 

2.2 Estimates on Characteristics 

In this section we investigate the linear Vlasov problem 

dtf + dwM"^} + div^/^a} + ^div^r (a - (3\v\^)v} = 0, (t, x, v)£R+xR^x R'^ (11) 

no) = r (12) 
where a E L~(M+; W^'°°{R'^)) and G 7W+(M^ x M"^). 

Definition 2.1 Assume that a e L°°(R+;W'^^'^{R'^)) and /'"^ e M'^{R'^ x R'^). We say 
that f^ G L°°(M+; A^f,(M^ x R'^)) is a measure solution of (fTT]) - (fT2|) if for any test function 
(f e C^{R+ xR'^ X M^) we have 

I I {dt + vV^ + a-V^ + -{a-(3\v\'^)vV^,}^f{t,x,v)d{x,v)dt 

+ / ip{0,x,v)f'''{x,v)d{x,v) =0. 

We introduce the characteristics of the field v ■ Vx + a ■ Vv + l{c( — P\v\'^)v ■ 

d dV^ 1 

~d^ = ^'(^)' ^ = «(^' ^'(^)) + -(« - /3 \V^s)\^)V%s) 



X^s = 0) =x, V^s = 0) =v. 

We will prove that {X"^ , V^) are well defined for any {s,x,v) G x M'^ x M'^. Indeed, on any 
interval [0,T] on which (X^ jV^) is well defined we get a bound 

sup {\X'{s)\ + \V'{s)\} < +00 

6G[0,T] 

implying that the characteristics are global in positive time. For that we write 

= ai^^X^s)) • V^s) + ^(a - l3\V'{s)\')\y%s)\'. (13) 

and then, we get the differential inequality 

< 2\\a\\L^\V%s)\ + -{a - P\V^s)f)\V^s)\^ 
as e 

for all s G [0,T], so that 

sup |y(s)| < +00, sup < l^l + T sup |y^(s)| < +00. 

se[o,r] se[o,T] sg[o,t] 

Once constructed the characteristics, it is easily seen how to obtain a measure solution for 
the Vlasov problem (jlip -(jl2p. It reduces to push forward the initial measure along the 
characteristics, see |10] for instance. 



Proposition 2.2 For any t G we denote by f^{t) the measure given by 

i;{x, v)r{t, X, v) d(x, v)= [ i;{{X',V'){t; 0, x, v))r{x, v) d(x, v) , (14) 



for all i/j G C°(IR'^ x R'^). Then the application t f%t), denoted /'"#(X^ y=)(t; 0, •, •) is 
the unique measure solution of (dH), 1^, belongs to C{R+;MbO^'^ x ^'^)) satisfies 

fit, x, v) d(x, v)= [ f "(x, v) d(x, v),teR+. 

Proof. The arguments are straightforward and are left to the reader. We only justify 
that f G C{R+;Mb{R'^ x M'^)) meaning that for any V G C°(IR'^ x R'^) the apphcation 
t —J" f^d-x^d'4'{x,v)f^{t,x,v) d{x,v) is continuous. Choose ip G C)?(M'^ x M^). Then, for any 
< ti < ^2 we have 



V'(x, v)f{t2, X, v) d(x, v) - tpix, v)f{ti,x, v) d{x, v) 

mX', V'){t2-M,x, v)) - i,{x, v)] f{t,,x, v) d(x, v). 
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Taking into account that {X^, V^) are locally bounded (in time, position, velocity) it is easily 
seen that for any compact set K C x M'^ there is a constant C{K) such that 

\X^{t2]ti,x, v)-x\ + \V'{t2; ti,x, v)-v\< \t2 - ti\C{K), (x, v) G K. 

Our conclusion follows easily using the uniform continuity of ip and that ||/^(ii)||A46 = 
||/™||a46- Notice also that the equality (fH|) holds true for any bounded continuous func- 
tion ip. ^ 

We intend to study the behavior of {f^)e when e becomes small. This will require a more 
detailed analysis of the characteristic flows {X'^,V^). The behavior of these characteristics 
depends on the roots of functions like A + ^{a — (3p'^)p, with p G M+, ^ G M. 

Proposition 2.3 Assume that A < and < e < 2ar/(|^|3\/3)- Then the equation 
A^(p) := eA + (a — f3p'^)p = has two zeros on M^, denoted pf{A),p2{A), satisfying 

0<pl<^<p2<r 



and 



lim^ = ^, lim!:^ = M 

e\0 £ a e\0 £ 2a 



where r = y^a/ f3. 



Proof. It is easily seen that the function A*" increases on [0,r/-v/3] and decreases on 
[r/^/S, +00 [ with change of sign on [0,r/^/3] and [r/^/3,r]. We can prove that (pf )e, (/02)e 
monotone with respect to e > 0. Take < £ < £ < 2ar /{\A\3^/3) and observe that A"^ > A^. 
In particular we have 

X%pl)< X^pl) = = x%pI) 
implying pf < pf, since A^ is strictly increasing on [0,r/\/3]. Similarly we have 

A^(p|) < A^(p|) = < X%pI) 



and thus pi > P21 since A^ is strictly decreasing on [r/\/3, r]. Passing to the limit in A^(p' 
0,k £ {1,2} it follows easily that 

lim pf = 0, lim po = r. 

Moreover we can write 

dp ^ e\o pf £\0 pf 



and 



dp e\o pI ~ ^ '^^'^ PI ~ 

saying that 

□ 

The case A > can be treated is a similar way and we obtain 

Proposition 2.4 Assume that ^ > and e > 0. Then the equation A^(p) := eA + (a — 
Pp'^)p = /ias one zero on denoted p^{A), satisfying 

^ Pl-r \A\ 
p%>r, \\m.^ ='—^. 

Using the sign of the function p — )• e||a||ioo + (a — /3p'^)p we obtain the following bound for 
the kinetic energy. 

Proposition 2.5 Assume that a G L°°(R+;W'^'°°(U.'^)), (1 + li-Hf" € x R'^) and 

/ei us denote by the unique measure solution of (jlip . ()12p . T/ien we have 



v\^f''{-,x,v)d{x,v) 



< I Mr + \v\v'''ix,v)d{x,v). 

-) 



Proof. We know that 



^AV'f < 2||a|U^|F^(t)| + -{a - P\V%t)\^)\V%t)\' = -\V%t)\X'{\V%t)\), t G M+. 
at e e 

By comparison with the solutions of the autonomous differential equation associated to the 
righthand side, we easily deduce that 

\V%t;0,x,v)\ < max{\v\,pl{\\a\\L-=)} , 

for any T G M+, {x,v) G x M'^. This yields the following bound for the kinetic energy 

\vfr{T, X, v) d{x, v)= [ \V'{T; 0, x, f )| v) d{x, v) 

„e\2 , L,|2i^in 



< / Mr + \v\v'''ix,v)dix,v). 

□ 

The object of the next result is to establish the stability of around \v\ = r. We will 
show that the characteristics starting at points with velocities inside an annulus of length 
proportional to e around the sphere rS get trapped there for all positive times for small e. 
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Proposition 2.6 Assume that e\\a\\Loc < 2ar/(3\/3) and that p2{—\\a\\Lao) < \v\ < /03(||a||L°°)- 
Then, for any {t, x) G M+ x M'^ we have 

pIHMl^) < \V%t;0,x,v)\ < pU\\a\\L^). 

Proof. As in previous proof, we know that 

±\V^\^<l\V^t)\X^\V%t)\), teM+. 

By comparison with the constant solution to the autonomous differential equation as- 
sociated to the righthand side, we get that sup^gjj^ \V'^{t;0,x,v)\ < p-^. Assume now that 
there is T > such that |y^(r)| < pi and we are done if we find a contradiction. Since 
|y^(0)| = 1^1 > P2, we can assume that minjg[o,T] 1^^(01 > Pi > by time continuity. Take 
now t G [0,T] a minimum point of t — t- on [0,T]. Obviously t > since 

\V'{t)\<\V%T)\<p^2<\v\ = \V'm- 

By estimating from below in (|13p and using that t is a minimum point of t — )• > on 

[0, T], we obtain 

..d {a - f3\V%m\V%t)\ X^lV^m 

U > —\V [t)\ > —\\a\\L^ H = . 

at e £ 

But the function has negative sign on [0, Pi]U[/32, +oo[. Since we know that min^gfg.T] 1^^(01 > 

pf, it remains that 

min \V%t)\ = \V'(t)\ > p% 
ie[o,T] ' ^ ^' ' ^ ^1 - 

which contradicts the assumption |y^(T)| < pi- □ 

Let us see now what happens when the initial velocity is outside [p|(— ||fl||L°°)) P3(||a||L°°)]- 
In particular we prove that if initially v ^ 0, then V^{t),t G M4. remains away from 0. We 
actually show that the characteristics starting away from zero speed but inside the sphere rS 
will increase their speed with respect to its initial value while those starting with a speed out- 
side the sphere rS will decrease their speed with respect to its initial value, all for sufficiently 
small e. 

Proposition 2.7 Consider e > such that e\\a\\L°° < 2ar/{3V3). 

1. Assume that < l^l < p|(~||fl||L°°)- Then for any {t,x) G x M'^ we have 

< \v\ < \V'{t;0,x,v)\ < p|(||a||Loc). 

2. Assume that /03(||a||L°o) < \v\. Then for any {t,x) G M*^ x M.'^ we have 

PK-lklli-) < \V%t;0,x,v)\ < \v\. 
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Proof. 1. Notice that if |y^(T; 0, x,v)\ = P2 for some T > 0, then we deduce by Proposition 
EJlthat /9| < \V%t)\ < pI for any t > T and thus \V^ {t]d,x,v)\ > pi > > T. It remains 
to estabhsh our statement for intervals [0,T] such that < pi fo^ ^'^Y ^ ^ [0;^]- We 

are done if we prove that t — t- is strictly increasing on [0,T]. For any r g]0,T] let us 

denote by t a maximum point of t — )• > on [0, r]. If t G [0, t[ we have < 

ctnd tlixis 

..d , {a-m^Cm\V^(i)\ x%\v^m 

> —\v (t)\ > — a 2,00 H = . 

at e £ 

By construction < pi ^-iid moreover, 

\V'(t)\ =max|y^| > \v\ > pI, 

[0,r] 

and thus, A^(|y^(t)|) > for all t £ [0,T]. Consequently, we infer that t — ?■ is strictly 

increasing on [0, T] since 

^lT/^mi> WnW , {a - f^\V%m\V%t)\ A^CZlM^n 

-r- r v*j ^ — a L°° H = > U. 

at £ e 

Therefore we have t = t saying that |T^^(t)| > \v\ for any r G [0,T]. 

2. As before, it is sufficient to work on intervals [0, T] such that > P3(||ct||L°°) for 

any t £ [0,T]. We are done if we prove that t — )• is strictly decreasing on [0,r]. We 

have for any t £ [0, T] 

d |....,.| . II II ^ {a-m%m\VHt)\ X'{\V%t)\) 
at £ £ 



where for the last inequality we have used > pl,t £ [0,T]. 



□ 



3 The limit model 

We investigate now the stability of the family {f^)s when £ becomes small. After extrac- 
tion of a sequence {£k)k converging to we can assume that converges weakly -k in 
L°°(M+; A^f,(R'^ X M"')), meaning that 

lim / / ip{t,x,v)f'^''{t,x,v)d{x,v) dt = / / ip{t,x,v)f{t,x,v)d{x,v)dt 

for any ip £ L^{R^; (M*^ x M.'^)). Using the weak formulation of (fTT]) - (fT2]) with test functions 
r]{t)ip{x,v), r] £ C<}(M+), £ C^iR'^ x R'^) one gets 

[ [ {r]'{t)ip + 7]{t)vV^ip + 7]{t)a-V^p}f>={t,x,v)dix,v)dt 

+— I I ri{t){a-f3\v\'^)v-Vy(pf''{t,x,v)d{x,v)dt 
Sfc jr+Jr^xM'' 

= -/ r]{0)(p{x,v)f''{x,v)d{x,v). 
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Multiplying by Ek and passing to the limit for k — >• +00 yields 

'n{t){a - I3\v\'^)v 

■ ^vff{t, X, v) d{x, v) dt = 
and therefore one gets for any t G M+ and ip G C^(M'^ x R'') 

(a — /3|?;p)v • Vvff{t, X, v) d(x, v) = 0. 



Under the hypothesis (1 + Iv]"^)/"^ G x M.'^) we deduce by Proposition [23] that (1 + 

^ A^^(M'^ X M'^) and therefore, applying the {x,v) version of Proposition [LT] (whose 
proof is detailed in the sequel) , we obtain 

supp f{t) CR'^ X ({0} U rS), t G ]R+. 

The proof of Proposition 11.11 is based on the resolution of the adjoint problem 

-{a- /3\v\'^)vV^ip = ^{v), vGM.'^ 

for any smooth righthand side ip with compact support in "^({0} U rS). 

Proof, (of Proposition fTTT]) It is easily seen that for any F G A^^(E'^x]R'^), supp F C {OjUrS 
and any ip G C^(M'^) we have 

(a - ^\v\^)v ■ V^ip{v)F{v) dv = 



saying that div„{F(a — /3|up)u} = 0. Assume now that div„{F(a — /3|f = for some 
F G M'liR''- X M'^) and let us prove that supp F C {0}UrS. We introduce the flow V = V(s; v) 
given by 

^ = {a-(3\V{s;v)f)V{s;v), V{0;v) = v. (15) 
ds 

A direct computation shows that |^ are left invariant 

(a - m')v ■ V. (^) = (a - /3|.P) (/ - ^) ^ = 

and therefore 

V{s;v) = \V{s;v)\^, v ^ 0. 

\v\ 

Multiplying ([15]) by V{s;v)/\V{s;v)\ yields 

^jV\ = {a-ms;v)\'ms;v)\ 



whose solution is given by 

\Vis;v)\ = \v\ 
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Finally one gets 



with S{v) = — oo if < |f I < r and S{v) = 2^ In ^1 — < if \v\ > r. Notice that the 
characteristics V{-]v) are well defined on M_|_ for any v G and we have 

V 

lim V(s; f ) = r— if f 7^ 0, lim V{s;v) = if = 

and 

lim |V(s)| = if < If I < r, lim |V(s)| = r if \v\ = r, lim |V(s)| = +cxd if \v\ > r. 
s\S{v) s\S(v) s\S{v) 

Let us consider a function ip = ip{v) with compact support in '^({0} U r§). We intend to 
construct a bounded function (p = ip{v) such that 

-{a- l3\v\'^)vVyip = i^{v), vGM'^. 

Obviously, if such a function exists, we may assume that v^(0) = 0. Motivated by the equality 

-^W{^{s]v))} =il)iV{s]v)), 0<|r;|<r, -00 < s < 
as 

and since we know that lims_i._oo V(s; f ) = for any < < r, we define 

ip{v) = -j V(V(r;u))dT, < |u| < r. (16) 

Let us check that the function ip in (jl6p is well defined and is in |?;| < r. The key point 
is that ij) has compact support in "^({0} U rS) and therefore there are < ri < r2 < r < < 
r4 < +00 such that supp <Z {v : ri < \v\ < U {f G R"^ : < \v\ < ri}. It is 

easily seen that r — t- |V(r;f)| is strictly increasing for any < |f| < r. Therefore, for any 
\v\ < ri we have |V(r;w)| < |V(0;f)| = |f| < ri, r < 0, implying that 

(p{v) = -[ V(V(r;v)) dr = 0, < \v\ < n. 



For any v with ri < < r2 there are ri < < T2 such that |V(ri;t))| = ri < r2 = |V(r2;t^) 
The time interval between ri and T2 comes easily by writing 



(a-/3|V(r)|2)|V(r)| 

implying that 

|''"2| + I'T'll = ^"2 — ■7"! = 



1 



'■2 dp 

J (a-/3p2)p' 
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From the equality 



?A(V(r; v)) At- V(V(r; v)) dr = - / ^/;(V(r; v)) dr . 

-OO J Tl J Tl 



we deduce that 

dp 



W{v)\ < Inl IIV^IIco < / ll^llco. (17) 



ri 



Assume now that r2 <\v\ < r. There is r2 > such that v = V(r2;r2|^) and therefore 

/■° f° V 

(p{v) = - I iI){V{t;v)) At = - I V'(V(r + T2; r2— r)) dr 



\v\ 



= -[ ipiVir + T2;r2^)) dr = - [ tp{V{T;r2^)) dr = ip ( r27^] . 

J -00 m J -00 m V mJ 

In particular, the restriction of (/? on r2 < |v| < r satisfies the same bound as in ()17p 

Mv)\< / 7 ^-^UWco, r2<\v\<r. 

It is easily seen that 93 is C"*^ on < |f | < r. For that it is sufficient to consider ri < < r2. 
Notice that 

kt;v) = — — 12 ) 



dv ■ 1^1 V ^ 

and therefore the gradient of ip remains bounded on ri < < r2 

V^^{v) = - J — (r; t.)V^(V(r; v)) dr 

since on the interval r G [ri, 0] we have | V(t; v)\ £ [ri, 1^1] C [^i, r2]. Taking now as definition 
for |t;| = r 

we obtain a bounded function on < r satisfying 

-(a - f5\v\^)v • V„<^ = Iv^WI < / ' 7 llV'lIco, \v\ < r. 

We proceed similarly in order to extend the above function for \v\ > r. We have for any s > 



—(p{V{s;v))+(p{v)= / iP{V{t;v)) dr, \v\ > r. 
Jo 

As lims^4_oo V(s; v) = r-^ we must take 

ip(v) = lim \ ip(V(s;v)) + / V(V(t; f )) dr > = ( r-^ ) + / ^(V(r;f))dr, l^;! > r. 
*-^+°° I Jo J V 1^1/ Jo 

Clearly, for any \v\ > r the function r — )• |V(r;t')| is strictly decreasing. Therefore, for any 
r < < r3 we have 

ifiv) = if (^r^^ = if (^r2^ 
15 



since |V(r;v)| < |f| < and V'(V(r;f)) = 0, r > 0. If r3 < |?;| < let us consider 
r4 < < T3 such that |V(r3;ti)| = < = |V(r4;f)|. The time interval between r4 and ra 
is given by 



ksl + \t4\ = Ts - r4 
and therefore one gets for < \v\ < r4 



J r4 



dp 



(a - f3p'^)p 



< +00 , 



l^(^^)l < 



< 



LP r 



r2 



+ 
dp 



^{V{t-v)) dr 



+ 



dp 



CO- 



n {a-f3p'^)p J^^ {a-/3p'^)p 
Consider now \v\ > r4^. There is r4 > such that ?^4^ = V{t4;v) implying that 



(18) 



ip{v) = ip[r—\ + 



4){V{t;v)) dr = If [ r— + 



Lf r 



i>{V{T;V{u;v)))dT = ip{r- 
V \v\ 



^(V(r;w)) dr 

^{V{T;n--)) dr 



<p r4 



We deduce that the restriction of ip on {v : \v\ > r4} satisfies the same bound as in ()18p . 
Moreover the function ip is on {v : \v\ > r}, with bounded derivatives. Indeed, it is 
sufficient to consider only the case < \v\ < r4, observing that 

VMv) = ^ (^- ^) V.V^ (-2^) + £ |^(r;.)VV'(V(r;.)) dr 



|V(r;u)| e [r3,\v\] C [r3,r4], r G [0,r3], jrs] + |r4| 



dp 



/r4 {a-f3p'^)p 

By construction we have —(a — • V^93 = ip{v), \v\ > r. 

Consider a C"^ decreasing function on such that xl[o,i] = 1) X[2,+oo 



< +00. 



0. We know that 



{a- P\v\^)vVAipiv)x 



R 



F{v)dv = 0, R>0, 



saying that 



X (^) (a - m')v ■ V,<^ F{v) dv + Ija - /3|^|2),^(^)M^ ' |^M^ f{v) dv = 0. 

Since ip and ip = —(a — /3|f p)u • V^ip are bounded and F has finite mass and kinetic energy, 
we can pass to the limit for R — )• +00, using the dominated convergence theorem. We obtain 
for any function tp, with compact support in '^({O} U rS) 

7p{v)F{v) dv = - {a- /3\v\'^)v ■ V^ip F{v) dv = 0. 



Actually the previous equality holds true for any continuous function ip with compact support 
in ''({0} U rS), since Jj^dF{v) dv < +00, so that supp F C {0} U rS. q 
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In order to obtain stability for {f'^'')k we need to avoid the unstable equilibrium f = 0. For 
that we assume that the initial support is away from zero speed: there is ro > (eventually 
small, let us say tq < r) such that 

supp f"" C {{x,v) G M'^ X M'^ : \v\ > ro}. (19) 

Proposition 3.1 Under the hypothesis (|19p we have for any e > small enough 

supp f{t) C {{x,v) eM.'^xR'^ : \v\ > ro}, t G M+. 

Proof. Take e > such that e||a||Loo < 2ar/(3\/3) and /Oi(— ||a||L°°) < ^q. For any continu- 
ous function ip = ijj{x,v) with compact support in M'^ x {t> : \v\ < tq} we have 

i>ix, v)fit, X, v) d(x, v)= [ i>{X%t; 0, X, v),V%t; 0, x, v))r{x, v) d{x, v) 

d 

\v\>ro 

yrix,v) dix,v). 

But for any \v\ > ro > pf we know by Proposition 12.71 that |I^^(t; 0, x,v)\ > \v\ > tq, implying 
that ip{X'^{t),V^{t)) = 0. Therefore one gets J^a~^^d'>Pix,v)f^{t,x,v)d{x,v) = saying that 
snpp f%t) C {{x,v) :\v\>ro}. □ 

We are ready now to establish the model satisfied by the limit measure /. The idea is to 
use the weak formulation of (jlip . ()12p with test functions which are constant along the flow 
of {a- f3\v\^)vV ^ , in order to get rid of the term in j . These functions are those depending 
on X and Surely, the invariants ^ have no continuous extensions in f = 0, but we will 
see that we can use it, since our measures vanish around v = 0. 

Proof, (of Theorem II .ip We already know that / satisfies (jT]). Actually, since supp /^(t) C 
{{x,v) : \v\ > ro},t G M-|-,e > 0, we deduce that supp fit) C {{x,v) : \v\ > tq} and finally 
supp f{t) c M'^ X rS, t G M+. We have to establish ([6]) and find the initial data. Consider a 
decreasing function x on M+ such that xl[o,i] = liX[2,+oo[ = 0- ^oi any r] = rj{t) G C^(M+), 
(p = ip{x,v) G C^{W^ X M.'^) we construct the test function 

'2\v\ 



9{t,X,v)=7]{t) 



1-X 



V 

(f I X, r — 

\v\ 



Notice that 9 is and 9 = for \v\ < When applying the weak formulation of (jlip - (fT2 
with 9, the term in - vanishes. Indeed, we can write 



rj{t){a- f3\v\^)v -V 



2\v\ 



\v\ 



V ( ] \ f%t, X, v) d(x, v) dt 



-[ ri{t) [ {a-(3\v\'^)vvAv(x,r^]} f{t,x,v)d{x,v)dt = 0. 
e Jr+ J\v\>ro I V \v\ 
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For the term containing dt9 we obtain the following limit when k — >■ +00 
r^:= / / dt9r>'{t,x,v)d{x,v)dt^ I I dt0fit,x,v)dix,v)dt 

= v'it) (p (x,r^] f{t,x,v)d{x,v) dt 

Jr+ J\v\>ro \ \V\J 

= v'{t) ^ix,r^] f{t,x,v)d{x,v) dt 

Jr+ J\v\=r \ m/ 

= / v'{t) / (p{x,v) f{t,x,v)d{x,v) dt 

Jr+ J\v\=r 

= / dtir](p)f{t,x,v)d{x,v)dt. 

jR+jRdxR'i 

Similarly, one gets 

:= [ [ vVx0f''it,x,v)dix,v)dt^[ [ v ■Vxef{t,x,v)d{x,v) dt 

jR+jR'^xR'i JR+jR'^xR'i 

= / / V ■Vx{rj(p)f{t,x,v)d{x,v) dt. 

jR+jR''xR'i 

For the term containing a • notice that on the set |v| > ro we have 

a ■ Vv9 = rj{t)a -Vvlipix, r-j-r > = r]{t)-r-;a ■ [I t-t^ (V^<y?) ( x, r- 



v\ / \ ' ~ ' l^'l V \v\^ / ' 'V \v 



and therefore we obtain 



T^:=f I a-VJf''{t,x,v)d{x,v)dt^ [ [ a ■ Vv9f{t,x,v) d{x,v) dt 

Jr+Jr''xR<'- Jr+Jr'^xR'' 



= [ vit) [ 

Jr+ J\u 



( / - ^^r^ ) a ■ {Vv(p) ( X, r-^ ) f{t, X, v) d{x, v) dt 
>ro I""! V \w ) \ H) 

I 1— ^ a ■ Vv{w)f{t, X, v) d{x, v) dt. 

\v\ / 



For treating the term involving the initial condition, we write 

T4 := [ 6(0, X, v)r{x, v) d{x, v)= [ r?(0)^ fx, r^) r{x, v) d{x, v) 

jR'lxR'l JR'lxR'i \ MJ 

= [ ri{OMx,v){r){x,v)dix,v). 

JR'^xR'i 

Passing to the limit for k — )• +00 in the weak formulation + T2 + T4 = yields the 

problem 

dtf + divMv} + div, |/ 1^1 - ^) a| = 0, /(O) = (/-) 



as desired. 



□ 
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Remark 3.1 The constraint ([7|) is propagated by the evolution equation This comes by 
the fact that the flow {X, V) associated to the field v ■ + ^/ — a ■ V„ leaves invariant 
R"^ X rS. Indeed, if{X,V) solves 



dX dV ( V{s)(^V{s) . 



X{s; 0, X, v) = X, V{s] 0,x,v) = v ^ 

then 



2ds' ' " V 

saying that \V{s;0,x,v)\ = \v\ for any {s,x) £ M-j. x W^. In particular, for any continuous 
function tp = ip{x, v) with compact support in "{W^ X rS) we have 



ilj{x,v)f{s,x,v)d{x,v) = / 'iIj{X{s;0,x,v),V{s;0,x,v)) (^f™^ {x,v)d{x,v) 

i;iX{s; 0, X, v), Vis; 0, x, v)) (f (x, v) d{x, v) = 



' |ti|=r 

since supp (/™) C M*^ x rS. Therefore for any s G M+ we have supp f{s) C M'^ x rS implying 
that diYy{f{s){a- p\v\'^)v} = 0,s G M+. 

Remark 3.2 By the uniqueness of the solution for ([6]) with initial data (/™), we deduce 
that all the family (/^)e converges weakly ★ in L°°(R+; Mb(R'^ x M*^)). 

4 The non linear problem 

Up to now we considered the stability of the hnear problems ()ll|) -(fT2|) for a given smooth 
field a = a{t, x) G -L°°(R_|_; W^'°°(R'^)). We concentrate now on the non linear problem 

dtf + div,{r^;} + div^ra^} + ^{fia - = 0, {t, x, v)eR+xR''x M'^ (20) 

with = —VxU-kp^ — H -k . The well posedness of the non linear equation (|2Up comes by 
fixed point arguments in suitable spaces of measures, and it has been discussed in |101 U] in 
the measure solution framework. We summarize next the properties of the solutions (/'^)e>o- 

Proposition 4.1 Assume h G Cl{W^),U G Cl{W^) and (1 + \v\'^)f''^ G 7W+(M^ x M^). For 
alle > 0, there is a unique solution (/^a^) G C(M+;Pi(M'^ x W^)) x L°°(M+; to 

dtf' + div^{/^t;} + dw^if'a'} + ^div4/"(a - ^\v\^)v} = 0, (t, x, G M+ x M"^ x M"^ (21) 



-VxUi^j fdv-Hi.f, H{x,v) = h{x)v (22) 
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with initial data /^(O) = satisfying the uniform bounds 



sup 

£>o,tm+ 



/ d(x,i;) < +00 



sup ||a^||L°°(R+;L^(Rd)) =: A < +00, sup \\'^xa^\\L°°(R+;L°°iR'i)) ='■ ^1 < 

£>Q £>0 

Moreover, if the initial condition satisfies 

supp C {{x,v) eM.'^ xW^ : |x| < Lo,ro < \v\ < Rq} 

for some Lq > 0, < tq < r < i?o < +00, then for any e > small enough we have 

supp/^(t) C G M'^ X M"* : \x\ < Lo + tRo,ro < \v\ < Rq}, t G M+. 

Proof. Here, we only justify the uniform bounds in e, the rest is a direct apphcation of the 
results in \10\ 0] . The divergence form of ()2ip guarantees the mass conservation 

[ fit, X, v) d(x, v)= [ f^{x, v) d{x, v), te M+. 

Notice that the term —divy{f^H-kf^} balances the momentum 

/ vdiv^{fHi<f'}d{x,v)= [ h{x-x'){v' -v)f{t,x',v')f{t,x,v) dix',v')d{x,v) = 

and decreases the kinetic energy 

v\'^divy{f^Hi.f'} d{x,v) = 2 [ h{x - x'){v' - v) ■ vf'{t,x',v')f%t,x,v) d{x' ,v')d{x,v) 

= - / h{x- x')\v - v'l'^fit, x', v')f{t, x, v) d{x', v')d{x, v). 

In particular, as G M'^(M.'^ x M'^), then the kinetic energy f^2d Iv]"^ f'^{t,x,v) d{x,v) 

remains bounded, uniformly in time t £ and e > 0. Indeed, using the continuity equation 
one gets 

V . {V,U ^ p')r{t,x,v) dix,v) = ~ [ iU*p'mx)p'{t,x)dx 
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and after multiplying (f2T|) by together with ([22]) . we obtain 

d r nvl^ Ujr£_\ ^^^^^ ^^^^ t;) - - / ia\v\'' - /3\v\')rit, x, v) d(x, v) 
dt Jmd^Tssd V 2 2 J e ' 



'xll 



h{x - x')\v - v'\^f%t, x', v')f%t, X, v) d{x', v')d{x, v) < . (23) 



Consider now t^ a maximum point on [0,T],T > 0, of the total energy 

W%t)= [ (\:^ + ^L!L^)f^(t,x,v)dix,v), tG[0,T]. 
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If = then it is easily seen that for any t £ [0, T] 

\^f{t,x,v) d{x,v) < [ ^-^r{x,v) d{x,v) + \\U\\l^ ( [ r ^{x,v)^ ' 
If g]0, r] then ^W^{t'^) > implying from ()23p by moment interpolation in that 



sup / {1 + \v\'^)f'^{f ,x,v) d{x,v) < +CO 

>0,T>0 JRdxRd 



e>o,: 

and thus the inequality Wit) <W%t^),t£ [0,T] yields 



sup / —f{t,x,v)d{x,v)< sup / — d(x,f) 

e>0,te[0,T] JR'^xRd £>0,T>0 JRdxR<* ^ 

+ f°d(x,7;)) <+oo. 

\JRdxRd / 

Therefore the kinetic energy remains bounded on [0,T], uniformly with respect to e > 0, 
and the bound does not depend on T > 0. The uniform bounds for come immediately by 
convolution with and H, thanks to the uniform estimate 

sup / \v\f^(t,x,v) < +00. 
£>0,teR+ JRdxRrf 

We analyze the support of (/^)e>o- Take e > small enough such that eA < 2ar/(3\/3) and 
pI{—A) < vq, p^{A) < Rq. By Proposition 13. II we already know that 

supp f{t) C {{x,v) eM.'^xR'^ : \v\ > tq}, t G M+. 

For any continuous function = il){x,v) with compact support in M'^ x {f G M*^ : > i?o} 
we have 



ij{x, v)f{t, X, v) d{x, v) = I ij{X'{t), V'{t))nx, v) d(x, v) 

^{X'{t), V'{t))l^ro<\v\<R,}r{x, V) d{x, V). 



We distinguish several cases: 

1. If ro < |^;| < P2i-A) we deduce by Proposition [22] that \v\ < \V^{t;0,x,v)\ < pKA) < 
Ro, t G M+,e > 0. 

2. If P2i-A) < \v\ < pI{A) we obtain by Proposition [Ml that pll"^) < 0, x, f )| < 
<Ro, te R+,e> 0. 

3. If pI(^) < \v\ < Ro one gets thanks to Proposition O pIC"^) < 0, x, w)| < 
\v\ < Ro. 

In all cases {X^ , V^){t; 0, x, v) remains outside the support of ip, implying that 

^p{x,v)f^{t,x,v) d{x,v) = 0. 
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Thus for any t G M+ and £ > small enough one gets 

supp f{t) C {{x,v) G X M'^ : tq < \v\ < Rq}. 

Consider 6 £ C^(M) non decreasing, verifying 6{u) = if u < 0, 6{u) > if n > 0. Applying 
the weak formulation of (|2ip - (j22p with the test function 0{\x\ — Lq — tRo) yields 



ei\x\ - Lo-tRo)rit, X, v) d{x, v) = / ei\x\ - Lo)rix, v) d(x, v) 

+ 11 e'{\x\-Lo-sR^))(v^-R^ f{s,x,v)<i{x,v)ds<Q 

Jo JR^'xRd. V fI / 



implying that supp f%t) C {{x, v) GR'^xW : \x\ <Lq + tRo},t £ 



□ 



The uniform bound for the total mass allows us to extract a sequence {£k)k C con- 
vergent to such that {f^^)k converges weakly * in L°°(M+; A^;,(M'^ x W^)). The treatment 
of the non linear term requires a little bit more, that is convergence in C(M-)- ; Pi (M"^ x W^)) 
or at least in C{[5, +oo[;Vi(W^ x W^)) for any 5 > Q. The key argument for establishing that 
is emphasized by the lemma 

Lemma 4.1 Consider e > small enough. 

1. For any {x,v) € M'^ x M*^ with ro < \v\ < P2{—A) — £, the first time tf = such that 
\V^{ti;0, x,v)\ = P2{—A) — £ satisfies 

t\ < ——IT In 

2. For any {x,v) ^W^xW^ with pK^) + £ < \v\ < Rq, the first time t^ = ^2(^''^) suc/i that 
|y^(t2; 0, X, u)| = p'^{A) + £ satisfies 

Proof. 1. During the time [0,tf] the velocity modulus remains in [tq, P2{—A) — e] C 

[pI{— A) , p2{— A)] and we can write for any t £ [0,tf] 



£ 



d\V^\ d\V^ 
dt ^ dt 



> 7^^^ = 1 



-eA+{a-l3\V^{t)\^) \V^{t)\ " a^^t,X^{t)) ■ + l{a - P\V^{t)\^) \V^{t)\ 

since —£A + (a — [iu^)u is positive for u € [p\{—A),p2{~A)]- Integrating with respect to 
t G [0, tf ] yields 

t\[x,v)<£\ -. ^ < £ -. ^ „, ■ 

-£A + {a - f3u^)u -£A + {a-^u^)u 
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Recall that P2{—A) is one of the roots of u — >• —eA + (a — j3v?')u and therefore a direct 
computation lead to 

-£A+{a-(3u^)u = P{p%-u)[u^ +upl+{plf -r^] > 2/3ro(p|-ti), u G [ro,/o|], e smah enough 
implying that 

2. During the time [0, tl] the velocity modulus remains in [p^{A)+e, Rq] C [pKA), +oo[ 

and we can write for any t £ [0, tl] 

^ dt > dt _ -j^ 

eA + {a-f3\V^{t)\^) \V^{t)\ " a^^t, X^{t)) ■ + l{a - (3\V^{t)\^) \V^{t)\ 

since eA + {a — (3u'^)u is negative for u G [/jK^), +oo[. Integrating with respect to t G [0,^1] 
yields 

Pl{A)+e rpl{A)+e 



' /|^| eA + {a - j3v?')u eA + {a - I5v?)u' 



tUx,v)<e I —. : ^^<e 

J\v\ 

By direct computation we obtain 

eA+{a-/3u^)u = -(3{u- pDlu^ +upl + {pl)'^ -r"^] < -2^r^(n-p|), u > pi, e small enough 
implying that 

p. , e du e ^ fRo-p%\ e , f Rq - r\ 

□ 

We intend to apply Arzela-Ascoli theorem in C(M+; ■Pi(M'^ x W^)) in order to extract a conver- 
gent sequence {f^'')k with limfc_j.+oo ^fc = 0. We need to establish the uniform equicontinuity 
of the family (/^)e>o- The argument below is essentially similar to arguments in |10] . 

Proposition 4.2 1. If the initial data is well prepared i.e., supp f° C {{x,v) G M'^ X M'^' : 
< -Z>0; I'wl = ^} then there is a constant C (not depending on t £ IR+,e > 0) such that 

Wi{f{t),f{s))<C\t-s\, t,seR+,e>0. 

2. //supp /™ C {ix,v) G M'^ X M'^ : |x| < LQ,rQ < \v\ < Rq} then there is a constant C (not 
depending on t & M+, e > 0) such that for any 6 > we can find es satisfying 

Wi{f^t),f{s))<C\t-s\, t,s>6, 0<e<es. 
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Proof. 1. Consider ip = ip{x,v) a Lipschitz function on x with Lip{ip) < 1. For any 
t,s £ M+, e > we have 

V{nt)-r{s))d{x,v) = [ {^{X^{t),V%t))-^{X^s),V^s))}r{x,v)d{x,v) 

< [ {\X^t) - X^{s)\ + \V%t) - y^(s)|}l||,|=,}rd(x,t;). 

Thanks to Proposition 12.61 we have for any {t,x,v) G x M'^ x rS 

-r ^ \V%t;0,x,v)\ -r ^ -r 
e ~ e ~ £ 

and it is easily seen, integrating the system of characteristics between s and t, that 



|X^(t;0,x,v) -X^(s;0,x,'y)| 



y^(r;0,x,v) dr 



< Ro\t - s\ 



and 

\V^{t]Q,x,v) - V%s]Q,x,v)\ < 

<\t- s\^A + I3{r + Ro)Ro max 
Our conclusion comes immediately by Propositions 12.31 E 



2. Consider 5 > and £5 small enough such that In {^—p^) < S, In {^^^) < S for 
< £ < £s- For any Lipschitz function if with Lip(99) < 1 and any t, s > 5 we have 



(^(r(t)-r(5))d(x,7;) 

For any (r, x) G . 



< / {\X^{t)-X^is)\+\V%t)-V%s)\}l{ro<\v\<Ro}r dix,v). 
X M'*, - e < |u| < /9|(A) + e we have by Propositions [213 EH 

/^K-A) - e < |y^(T;0,x,^;)| < p^A) + £. 



The same conclusion holds true for any r > 5, x G M"* and \v\ G [ro, P2(~^)~£[U]/03(^)+e, -Ro]) 
thanks to Lemma HTTl since 5 > max{tl{x,v),t2{x,v)} (after a time 5, the velocity modulus 
|y^(r; 0, X, v)| is already in the set {w : P2{—A) — £ < \w\ < p^{A) +e}). Our statement 
follows as before, integrating the system of characteristics between s and t. \ — | 

Applying Arzela-Ascoli theorem, we deduce that there is a sequence {£k)k C M^, convergent 
to such that 



lim Wi{f'^{t),f{t)) = uniformly for t G [0, T], T > 



for some / G 



ViiR'^ X M'^')) if supp f° C {{x,v) G M"^ X M'^ : |x| < Lo,\v\ = r} and 



lim Wiif^ {t)J{t)) = uniformly for t G [5, T], T>6>0 
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for some 

It is easily seen that if the initial condition is well prepared then there is a constant C cf. 
Proposition 14.21 such that Wi{f{t),f{s)) < C\t — s\, t,s £ R+. The same is true for not 
prepared initial conditions /™. Take 6 > and es as in Proposition 14.21 For any < e < es 
we have Wi{f^{t), f^{s)) < C\t — s|, t,s > 6. For k large enough we have Sk < £s and 
therefore /^''(s)) < C\t — s\, t,s > 5. Passing to the limit as k goes to infinity 

yields /(s)) < C\t — s|, t,s > 5. Since the constant C does not depend on 5 one 

gets 

Wi{f{t),f{s))<C\t-s\, t,s>0. 

In particular we deduce that / has a limit as t goes to since {Vi{R'^ x R'^), Wi) is a complete 
metric space and therefore we can extend / by continuity at t = 0. The extended function, 
still denoted by /, belongs to C(R+;Vi{R'^ xR'^)) and satisfies 

Wi{f{t),f{s))<C\t-s\, t,seR+. 

The above convergence allows us to handle the non linear terms. We use the following 
standard argument [24:\ [TO] . 

Lemma 4.2 Consider f,g£ Vi{R'^ x R'^) compactly supported supp / U supp g C {{x,v) £ 
R'^ xR'^ : \x\ < L, \v\ < R}, and let us consider 

af = -VxU-k fdv-H-kf, ag = -VxU-k gdv-H-kg. 

Then we have 

\\af - agh^^^s^Bn) < {W'^lUh^ + (ll^lli- + ^R^VMloof^} Wiif,g) 
where Bji stands for the closed ball in R'^ of center and radius R. 

Proof. Take vr to be a optimal transportation plan between / and g. Then for any x G M'^ 
we have, using the marginals of vr 

\{VxUkf){x)-{VxUkg){x)\= [ VxU{x-x'){fix',v')-g{x',v')}d{x',v') 

= / / [VxUix-x')-VxU{x-x")]dTr{x',v',x",v"] 

< ||V^?7||loo /" [ \x' - x"\dTT{x',v',x",v") 

JR'^xR'^/RdxRd 

< \\VlU\\L^Wi{f,g). 
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The estimate for H -k f — H -k g follows similarly observing that on the support of tt, which is 
included in {{x\v' ,x'\v") G M^°' : \v'\ < R, \v"\ < R} we have 

\h{x - x'){v - v') - h{x - x"){v - v")\ 

< \h{x - x'){v" - v')\ + \h{x - x') - h{x - x")\ \v - v"\ 

< {\\h\\l^ + AR'WVMl^f iW - x'f + \v' - v'ff . 

□ 

We are ready now to prove Theorem 11.21 

Proof, (of Theorem ll.2|) The arguments are the same as those in the proof of Theorem 11.11 
except for the treatment of the non linear terms. We only concentrate on it. Consider {f^'')k 
with limfc_^+ooefc = such that limfc^+oo /(O) = uniformly for t £ [0,T],T > 

if supp /™ C {{x,v) : |x| < LoJ^^I = ">"} and limi^^^^Wi{f'^''{t), f{t)) = uniformly for 
t G [5,T],T > 5 > if supp /™ C {{x,v) : |x| < LQ,rQ < \v\ < Rq} for some function 
/ G C(IR+ ; Pi (M'^ X M"^)). Thanks to Proposition 11.11 we deduce (for both prepared or not 
initial data) that 

supp f{t) C {(x, v) £R'^ xR"^ : \v\ =r}, t> 0. 

The previous statement holds also true at t = 0, by the continuity of /. The time evolution 
for the limit / comes by using the particular test functions 



9it,X,v)=7]{t) 



1 / 2IH 

1 -X — 



LD X, r — 
\v\ 



with 7] G Ci(M+), if G Cl( 



). From now on we consider only the not prepared initial 



data case (the other case is simpler). We recall the notation = —VxU * f^dv — H-kf^ 
and we introduce a = —VxU * J^a f dv — H -k f. Since / satisfies the same bounds as (/^)e, 
we deduce that ||a||Loo < A, \\Vxa\\L°° ^ ^i- For any 5 > we can write 



/ [ {a'^ - V.e -a-VJ f}d{x,v)dt 



< 



V^ef^ d{x,v)dt 



+ 



f'8 /" /'+00 f' 

I I a-V^e f d{x,v)dt + / / {a^^ -V^e f" - a-V^e f} d{x,v)dt 



<2A6\\V^e\ 



CO 



rdix,v) + 



(a' 



^v0 l{|^;|<Ro}/'* d{x,v)dt 



+ 



/ / a-V,eir^-f)dix,v)dt 



(24) 



We keep 5 > fixed and we pass to the limit when k goes to infinity. Lemma 14.21 implies 
that the second term in the last right hand side can be estimated as 

Ik^" -a|lL°°(R'*xSflo) = hr^ -afh^-iRdxEn^) < CiRoWi{f''{t),f{t)) when k +oo 
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uniformly for t £ [5, T] , implying, for T large enough 



+00 



(a^fc -a)-V.,ef'' d{x,v)dt 



< c{Ro)\\e\\ci rw,{f^{t)j{t)) dt^O 



\v\<Ro 

when k goes to infinity. For the third term in the right hand side of ()24p we use the weak -k 
convergence lim.k_^_^.^ f'^''{t) = f{t) in x M'^) for any t > 6, cf. Proposition 12.11 

lim [ a-VJ{f^{t)-f{t))d{x,v) = 0,t>6 

and we conclude by the Lebesgue dominated convergence theorem 

a-V^e{f''{t,x,v) - f{t,x,v)) d{x,v)dt = 0. 

Passing to the limit in (|24p when k goes to infinity, we obtain 



lim 

fc— >-+oo Jg 



[ [ {a'" ■ VvOf" - a ■ VvOf} d(x, v)dt 
Sending 5 to we obtain that 



lim sup 

fc— >+oo 



< 2^(5 1 1 V^, 6*110-0 



lim 



0^*= • Vy0 f' d(x, v)dt 



o- • Vt,0 / d(x, v)dt . 



□ 



5 Diffusion models 

We intend to introduce a formalism which will allow us to investigate in a simpler manner 
the asymptotic behavior of (jlip and Q. This method comes from gyrokinetic models in 
plasma physics: when studying the magnetic confinement we are looking for averaged models 
with respect to the fast motion of particles around the magnetic lines. The analysis relies 
on the notion of gyro-average operator [B], which is a projection onto the space of slow 
time depending functions. In other words, projecting means smoothing out the fluctuations 
with respect to the fast time variable, corresponding to the high cyclotronic frequency. This 
projection appears like a gyro-average operator. Here the arguments are developed at a 
formal level. 

We first introduce rigorously the projected measure on the sphere rS for general measures. 
Let / G X W^) be a non negative bounded measure on x M'^. We denote by (/) 

the measure corresponding to the linear application 

«..«)W(..«)d(.v.)./ *(.,.^)w(.,.)d(.,,.). 

JKdxRd JRdxRd V \V\ J 
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for all V e C°(M'^ x R'^), i.e., 

/ i^{x,v){f){x,v)d{x,v) = ip{x,v)f{x,v)d{x,v) + ( x,r^ \ f{x,v) d{x,v) , 

Jr'^x'S.'' Jv=0 Jv^O \ rl/ 

for all V' G C^(M'^ x M'*). Observe that (/) is a non negative bounded measure, 

if) {x,v)d{x,v) = f{x,v)d{x,v), 

with supp (/) C M'' X ({0} U r§). We have the following characterization. 

Proposition 5.1 Assume that f is a non negative bounded measure on x W^. Then (/) 
is the unique measure F satisfying supp F cW^ x ({0} U rS), 

/ 'ilj(x,r-^] F{x,v)d{x,v)= [ ( x,r-^] f{x,v) d{x,v), e C^{M.'^ x R'^) 
andF = f on x {0}. 

Proof. The measure (/) defined before satisfies the above characterization. Indeed, supp (/) C 
X ({0} U rS). Taking now il){x,v) = ip{x)x{\v\/5) with ip G C°(M'^) and 5 > one gets 



if) {x,v)d{x,v) = / ip(x)f{x,v)d{x,v) 



v=0 



Passing to the limit for (5 \ yields 



+ J ^ix)x(^-j^ fix,v)d{x,v). 



/ V>{x) if) {x,v)d{x,v) = ip{x)f{x,v)d{x,v), if e C[ 

Jv=0 Jv=0 

meaning that (/) = / on x {0}. Therefore one gets for any t/; G C^(R^ x 
/ ix,r-^] if) {x,v)d{x,v) = t/j{x,v) if) {x,v)d{x,v) 

Jv+Q V rl/ J\v\=r 

= / fpix,v) if) {x,v)d{x,v) 



tp if) (x, v) d{x, v) - ^ if) {x, v) d{x, v) 

<Rd Jv=0 

= '(P{f)ix,v)d{x,v)- il;f{x,v)d{x,v) 

JR'ixR'i Jv=0 

= ix,r^] f{x,v)d{x,v). 
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Conversely, let us check that the above characterization exactly defines the measure (/). For 
any V G C°(M°' x M"') we have 

/ 'iIj{x,v)F{x,v) d{x,v) = / il^F{x,v) d{x,v) + / 'il)F{x,v) d{x,v) 

Jm'^xR'* Jv=0 Jvj^O 

■ip{x,v)f{x,v)d{x,v)+ / ip ( x,r—\ F{x,v) d{x,v) 
V=0 Jvyto V M J 

^l^{x,v)f{x,v)d{x,v)+ / ip ( x,r—\ f{x,v) d{x,v) 
V=0 Jvyto \ M J 

saying that F = (/). , 

By Proposition 15.11 it is clear that (•) leaves invariant the measures with support in 
M'^ U ({0} U rS) . Consider f e {R'^ xR'^). We say that div^{/ (a - P\v\^)v} e MbiR'^ xR'^) 
if and only if there is a constant C > such that 



'xll 



{a- p\v\')vV^^f{x,v)d{x,v) < CIIV'IIloo, e x M"). 

In this case there is a bounded measure fi such that 

{a- ^\v\'^)vV^'ipfix,v)d{x,v) = [ ^pn, £ Cl{R'^ x R'^). 

By definition we take divt,{/(a — /3|f = fi. The main motivation for the construction of 
the projection (•) is the following result. 



Proposition 5.2 For any f G x R'^) such that div^{/(Q - P\v\'^)v} G Mb{R'^ x 

we have (divt,{/(a — = 0. 

Proof. Let us take divt,{/(a — /3|f p)^} = fi. We will check that the zero measure satisfies 
the characterization of (fi) in Proposition l5.ll Clearly supp = C M'^ x ({0}UrS). For any 
(fix) G C°(M'^) we have 

ip{x)fx(x,v) d{x,v) = lira f ip{x)x (^-^1 fJ-{x,v) d{x,v) 



lim I ^{x)x' ( ^ ) ^(a - P\v\^)f{x,v)d{x,v) = 



by dominated convergence, since 



< Q sup \x '{u)u\ + /3(5^ sup \x' {u)v? 



Therefore we deduce that div^{/(a—/3|t;p)t'} = OonM'^x{0}. Consider now V' G 
and lets us compute 

V' f X, r— ) f ) d(x, = lim / 1 3;, r— ) ( 1 — % ( -LJ- ) ) /i(x, v) d(x, 

v^O V 1^1/ '5\0yiRdxRd V 1^1/ V \^JJ 



lim 

5\0 



x,r— J x' y^j ~ l^\v?)f{x,v)d{x,v) = 
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since v • Vy{ilj{x, ^■^)} = 0. By density, the same conclusion holds true for any ijj £ 



X 



i"^) and thus (div„{/(a - /3\v\^)v}) =0. □ 



Remark 5.1 When f G A^^(M x M") does not charge M x {0}, (/) is given by 



or equivalently 



(/) = / V ( X, ) /, ^ G CO(M'^ X E'^). (25) 



Using Proposition 15.21 we can obtain, at least formally, the limit model satisfied by / = 
lime\^o/^- By Q we know that supp / C M*^ x ({0} U rS). The time evolution of / comes 
by eliminating f^^'> in ([5]). For that it is sufficient to project on the subspace of the measures 
satisfying the constraint (HD, i.e., to apply (•). 

{dtf) + (div,{/7;}) + (div,{/a}) = 0. (26) 

It is easily seen that {dtf) = dt (/) = dtf since supp / C M'^ x ({0} U rS) and therefore 
(/) = /. We need to compute the last two terms in (p6|) . We show that 

Proposition 5.3 Assume that a = a{x) is a hounded continuous field. Then we have the 
following equalities 

(div,{/^;}) = diY^ifv} if supp / C M'^ X ({0} U rS) 

(div4/a}) = div„ {/ - a} if supp / C E"^ x rS. 

As a consequence, ()26p yields the transport equation ([6|) obtained rigorously in Theorems 



Proof. For any tp G C^{W^ x M'^) we have 



V' (div^{/f }) = / ipdiv^ifv} + / i{j [x,r— ] div^ifv} 

id Jv=0 Jvj^O V 1^1 



lini [ tpx ( 4^ ) div^{/t;} + lini [ V (x,r-^) ( 1 - X ( 4^ ) ) div^{/i;} 
-lim/" • VajV-X ( ^ ) / - lim / v -V ( x,r^\ (l - x\^-t\] f 



V ■ ^xi^f = / il^diwxifv] 
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saying that {diYxifv}) = div^{/f }. Assume now that supp / C M'^ x rS. It is easily seen that 
divi,(/a) does not charge M'^ x {0}. Indeed, for any £ C^(M'^ x M.'^) we have by dominated 
convergence 



V'div^(/a) = lini / i'X [ -r ] div^(/«) 
v=0 "^NUjjjdxKd \ 



lim [ a ■ Vyipx (^-ri f ~ 1™ / ^ " 7~7T/t ' (^-71 ^/ = 0- 



Therefore we can use ([25 



-0 (div^(/a)) = / ^ [x, r— div„(/a) 
Iv^o V \v\ 



lim 

<5\o 



: — lim 

5\0 

+ hm 

<5\0 



1 - X ( -y j j -0 (^''^0 ) divt,(/a) 



--(f R -(-^^) -r)/ 



V' div, / / 



V (i9v 



□ 

We investigate now the limit when e \ of the diffusion model ([9]). We are done if we 
compute (At,/) for a non negative bounded measure with support contained in M*^ x rS. As 
before we can check that A^f does not charge M'^ x {0} and therefore, thanks to ()25p . we 
obtain after some computations 



i,{AJ)= V x,r— A,/= / A,^ x,r— /, VeC-, 



(27) 



Lemma 5.1 For any function (f € C (M \ {0}) and any r > we have 



Combining (|27p . Lemma l5. II and the fact that supp we obtain 



^ix,v){A,f) 



V V 



■.dli^{x,v)-2 



V ■ Vylp{x, V) 



ip{x, f )div^ < div„ 



f I 



V V 



+ 2/t 
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We deduce the formula 



(A^/) = div^ \ div^, 



f I 



V V 



+ 2/^ 



for any / satisfying supp / C M"^ x rS and the limit of the Vicsek model ([9]) when e \ 
becomes 



9t/ + div,(/z;)+div,<^/ / 



V (i9 V 



divt, < divt, 



f I 



V ^ V 



+ 2/i 



(28) 



with the initial condition /(O) = (/'"), as stated in ([10 



A Spherical coordinates and the Laplace-Beltrami operator 

In this appendix, we show the computations to relate the equations written in original vari- 
ables {x,v) to the equations in spherical coordinates {x,uj). Our limit densities have their 
support contained in M'* X rS and thus reduce to measures on R"^ x rS. For example, let us 
consider the measure on M'^ x rS still denoted by /, given by 



tp{x, u))f(x, uj) d{x, cj) 



for any function G C^(M'^xrS). In particular, to any / G A^^(M"'xM'*) not charging M'^x {0} 



Ip X, r— fix, v) dfx, v) 

v+Q \ \v\) 



it corresponds (/) G M.\ 



with supp (/) C M'^ X rS, whose characterization is 



^ x,r— f(x,v)d{x,v). 

v^O V \v\ 



We intend to write the previous limit models (in Theorems II. H 11.21 arid (j28p ) in spherical 
coordinates. 



Proposition A.l Assume that f G A^^(M x M*^), supp / C M x rS and let us denote by 



F G A^^(M^ X rS) its corresponding measure on M.'^ x rS. Therefore we have 



{dwMv)) = div,.(Fw), (div4/a)) = div^ |f (^/ - l(a; a [> , {A,f) = A^F 
Proof. Thanks to Proposition 15.31 we have for any ip G C^(M.'^ x rS) 



V'(x,w) {div^{fv)) 



Ip x,r— diva;(/f) 

v^O V \v\ 



V I V 

r--l ■ Vx'ip x,r— ] f 



V ■ V^V x,r— f 
vj^O V \v\ 



uj • Vxipix, uj)F 
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and thus {divxifv)) = diVa;(Fw). Similarly we can write 



/ — ^{u} (E)uj)] a ■ ( I — ^(w (g) cj) ) Vyip{x,uj)F 

'xrS\ ^ / V / 

/ - ^{u} (g) ojU a ■ \7^'ip{x,uj)F 
V'(x,a;)div^ I ( — (g w) ) a 



(div^(/a)) = div,^ <j F ( / - ^(o; w) ) a 



and therefore 



Here div^j stands for the divergence along rS (notice that (/ — ^{oj ^ w)) a is a tangent field 
of rS) and V^j = (yl — ^{u: ® w)) V^, is the gradient along rS. For the last assertion we appeal 
to the following well known result asserting that the Laplace-Beltrami operator coincides with 
the Laplacian of the degree zero homogeneous extension, see also [1]. 

Proposition A. 2 Consider ip = ip{uj) a function on r§ and we denote by ^ = ^{v) its 
degree zero homogeneous extension on M'^ \ {0} 

Therefore we have for any u £ rS 

Let us come back to the proof of Pr op osit ion I A . 1 [ For any ^ S C^(M'^ x rS) we introduce 



its degree zero homogeneous extension ^'(x,f) = ijj (x,rA). Thanks to Proposition IA.2I we 



can write 



/ ^{x,oj){A,f)= [ i,(x,r^){AJ)= [ ^{x,v)A,f=f A,^ f 

= / A^tP{x,v)f = / A^i;{x,uj)F= / i^{x,uj)A^F 

J\v\=r JR'^xrS Jr-^xtS 

meaning that (A^f) = A^F. |— I 
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For the sake of completeness, we finally write the equations in spherical coordinates in 
M'^. We introduce the spherical coordinates w = r(cos 9 cos 99, cos 9 sin ip, sin 9) with the angle 
variables {9, ip) g] — 7r/2, 7r/2[x [0, 2tt[, and the orthogonal basis of the tangent space to r§ 



eg = (— sin 9 cos if, — sin 9 sin 99, cos 9), e<^ = (— cos 9 sin 93, cos cos if, 0) 

with I eg I = 1, |e^| = cos 9. For any smooth function u on r§ we have 

el 1 

V^n = {VuiU • 69)60 + (V^u • 6^) — ^ = -deu eg H o^f^v^w 

cos^ 9 r r cos^ 



and for any smooth tangent field ^ = ^560 + ^,pe^ we have 

div^^ = - I ^-^^^(Ce cos 6*) + 5^^^ 
r [ cos 

The coordinates of the tangent field := F (/ — :pj{uj w)) a are ■Ce = C " = -Pogi 
-^^^ = Fa^o and we obtain 

cos"" ^ 









CO) 













-de{Fag cos 9) + d^{Fa^) 



The spherical Laplacian is given by 

r [cosy at' \ r / dip yrcos^y 
^ ^ ^ {^cos9deF) + -^BlF 



1 cos 6* 89 cos^ ''^ 



Proposition A. 3 The limit transport equation obtained in (|'28|) /or 



IS 



r [ COS J [ COS 9 o9 cos^ ^ 



We recall here the proof of Proposition IA.2I It is a consequence of a more general result. 

Proposition A. 4 Let us consider a function p = p{v) G C'^iW^), d > 2 which writes in 
polar coordinates p{v) = ip{p,a), p = \v\ > 0, o" = ^ G S. Therefore for any v ^ we have 

1 d 1 V 

^vV{v) = ^TZT^ip^'^dpip) + ^A„ip{p,a), p=\v\>0, cr = -r-r. 

p" '-op p^ \v\ 

Proof. Consider a smooth function ip = il){v) G with compact support in M^\{0}, which 
writes in polar coordinates ip{v) = tp{p,a), p = \v\ > 0, a = G §. We have 

— = Vvp ■ a, VvP = C^vf ■ cr)a + {I - a cr)Vv^ ^ 'dp ^ ~^ ^ ^=pa<^ 
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and 

dtp dip 

— = Vvip ■ <J, ^vW = iyv'4' ■ o-jo" + (/ - (g> a)Vv'p = + ^'^=p<tV- 

Integrating by parts yields 

Jm-i JRd Jm.+ [ op op P J 

= - / / 4 dp - / / « A,^ d. dp 

Js'^-^Jr+ dp\ op J J^^ p^ JsN-i 

and therefore 

19 1 V 

'^v^iv) = ^n^j7-{p^~'^dp'f) + ^A„(f{p,a), p=|'t;|>0, cr=-r-r- 
p" ^ op 

□ 

Proof, (of Proposition IA.2| ) The degree zero homogeneous extension <^{v) = Lp does 
not depend on the polar radius ^{v) = ^{(j) = f{uj = ra), a = j^. Thanks to Proposition 
\AA\ we deduce = -^Afj^ = ^A^^ip. Taking p = r, which means v = ra = uj we obtain 
A^,^>(a;) = A^ip{uj), oj £ r§. I 
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